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00 Abstract 

We suggest to use for XY2 molecules some results previously established in a series of articles for vibrational modes and 
electronic states with an E symmetry type. We first summarize the formalism for the standard u(2) D su(2) D so(2) 



in 



chain which, for its most part, can be kept for the study of both stretching and bending modes of XY% molecules. Next 
the also standard chain u(3) Z> u(2) D su(2) D so(2) which is necessary, within the considered approach, is introduced 
for the stretching modes. All operators acting within the irreducible representation (irrep) [NOO] = [NO] of u(3) are 
built and their matrix elements computed within the standard basis. All stretch-bend interaction operators taking into 

ft, 

g account the polyad structure associated with a resonance uji m L03 ~ 2 L02 are obtained. As an illustration, an application 
to the D2S molecular system is considered, especially the symmetrization in C^v It is shown that our unitary formalism 

43 ■ 

rj allows to reproduce in an extremely satisfactory way all the experimental data up to the dissociation limit. 
O 
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The hydrogen sulfide molecule and its isotopic species 
are of interest for terrestrial atmospheric pollutant mea- 
surements. As this gas is more heavy than air, it remains 
concentrated to the floor level and can be lethal at high 
concentration. Global warming process has increased the 
number of studies devoted to the chemical and physical 
properties of H2S and isotopic species [j], Q, S]. These 
molecules have also been observed in planet atmospheres 
like Jupiter or Venus JfSland appear in the analysis of 
the interstellar medium [□, il S] • 

Many papers, devoted to the analysis of the rovibrational 
spectra of D2S, have been published during the last thirty 
years : anharmonicity corrections to observed fundamental 
frequencies of vibration in 9] , various molecular structures 



have been evaluated in [l(J, [ll| for the V2 band, 
three fundamental bands 1^2 , v\ and v% of D2 32 S and the V2 
band of D 2 3i S, application of the MORBID (Morse Oscil- 
lator Rigid Render Internal Dynamics) computer program 
for the four isotopic molecules H 2 32 S, D 2 32 S , HD 32 S, 
and£T 2 34 S IjJ. 

Different theoretical models have been elaborated to im- 
prove the analysis of D2S or bent XY2 molecules in 



general 



R 



1 €31 ] . Initiated by the early works of Iachcllo 



and Oss [16(, algebraic formalisms seem to be good tools 
for the description of XYz molecules [ljj particularly to 
take into account the local behavior. Also the vibrational 
spectrum must be described using a Hamiltonian where 
the importance of Fermi-type interaction is taken into ac- 
count [18( . Other studies using potential energy surface 



have been recently published 



19 



20(. This last method is 
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particularly adapted when one has many data coming from 
different isotopic species. For further references about re- 
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cent analysis of H2S molecule and isotopes, the reader is 
invited to examine the references given in 

HQ, Ell- 
in this paper, we use and adapt formalisms that we devel- 
oped previously. Some of them dealt with unitary algebras 
applied in molecular spectroscopy 
studies 23| or works concerning al gebr aic chains which 



22j, pure algebraic 



24 



25 



may be adapted to XY 2 systems 
In the theoretical part we use the properties of two alge- 
braic chains, u(2) D su(2) D so(2) and m(3) D u(2) D 
su(2) D so(2), to analyze the bending and stretching vi- 
brational modes of bent XY2 systems. Basis states and 
operators oriented in these chains are built and the matrix 
elements of the oriented operators are given. Next sym- 
mctrization of these tensors in C*2u allows to build Hamil- 
tonian and tensor operators adapted to the C^t, molecular 
point group. A method to select, in the u(3) £g> u(2) dy- 
namical algebra, all relevant operators for a given polyad 
structure is proposed and applied to the case of a 2:1 res- 
onance. 

As an illustration, in the last section our approach is tested 
upon the D2S molecule. The experimental data are repro- 
duced with a standard deviation close to 0.5 cm -1 and the 
calculated dissociation energy is found close to the exper- 
imental one. 

2. Theoretical frame 

2.1. The standard u(2) D su(2) D so(2) chain 

We consider first physical systems under the assump- 
tion that they can be described from two boson creation 
and annihilation operators denoted bf and hi (i = 1,2) 
which satisfy the usual Bose commutation relations. A 
basis for the space of states may be obtained through re- 
peated action of b^ and b 2 onto the vacuum state |0, 0): 



K,n 2 ) = {7i 1 \n 2 \)- 1/2 bf ni ^" 2 |0,0>. 



(1) 



We note that these states may be taken as those of a two di- 
mensional oscillator (within a usual approach) or as those 



of two one dimensional oscillators. They will also be as- 
sociated with the dynamical states for a one dimensional 
oscillator. 

Tensor operators are built from the well-known Schwinger's 
realization of su(2) in terms of two boson operators [271 ]: 



J+ = b+b 2 , J- = b+h 



Jz = h{Ni-N 2 ) = \{btb 
with commutation relations 



1 -Kb), 



(2) 



[J„J±]=±J± , [J+,J-]=2J„ 

and u{2) is obtained with the addition of the linear invari- 
ant N = N ± +N 2 , with Ni |ni,n 2 ) = n t \ni,n 2 ) (i = 1,2). 

2.1.1. General tensor operators and states within the stan- 
dard chain 

Keeping with previous conventions 24| covariant su(2) 



states \jm)) and operators Tm are characterized by the 
relations (m : —j, ■ ■ ■ 



(3) 



Jz \jm)) = -to \ jm)) 

J± \jm)) =-[(i±m)(iTm + l)] 1 ' /2 IjtotI)), 
and likewise for irreducible tensor operators (ITO) 

\J zi Tm ] — TO T rn 

[J±, Tm] = -[(j ± m)(j T m + 1)]^ T^ v 

In the whole u{2) D su(2) D so(2) chain, symmetrized 
states and tensor operators are further characterized by an 
additional u(2) label [TO12 TO22] = [mi ~ TO2] m Gel'fand 



(4) 



notation 



28 



29] with j = (mi + to,2)/2 and 



'T 



Uh 



mi2 lmi-m 2 ] T U) j _ ^ g. 



(5) 



With equations ([2j |4j) it can be checked that a realization 
for the extremal components appearing in ([5]) is given by: 



T 



= (-l) m2 ! mi [ mi !m 2 !]-5 b+ mi b 2 m ' 2 (6) 



TOl — TO2 

mi 

This allows to generate two sets of fundamental operators 

[mi o ]r( i) = _ m)!( ^l + m)! ]-i/2 

(7) 



x 



[0- ro2 ] T (i) = ( _ 1)m2 [( ^2 _ -2 + m)!] _ 1/2 

2 2 (8) 

x b 1 ( ^ +rn ' > b 2 { ^- m \ 

with j — m2/2. In particular the standard covariant basis 
is obtained with ([7J acting upon the vacuum state: 

|[nO]jm)) = [n0] T^ |0,0 >= (-l)^ +m \j-m, j+m)). (9) 

^From the previous set ([71 [8]) one may build all functionally 
independent operators which may act within the irrep [n 0] 
of u(2) through 

r i U) 

[mi — m,2]q-(J) _ jj-jmax [mi 0] rp(-^) x [0 -m 2 ] J>{^) 

L Jm 

with 

|mi - m 2 | . toi + m 2 _ . 

Jmin — ^ _ J _ ^ Jmax ■ 

Their expansion in normal ordered form can easily be ob- 
tained; alternatively they may be written 

[mi -m 2 ] T (J) _ {mim 2 }g^ Ni+N2 j [m{ -m' 2 ] T (j) ^ ^ 

where 



/ _ mi-m 2 , ■ 
I I — 2 + J 



/ _ m 2 -mi , ■ 

h — ?. ' Ji 



hence j = (m' 1 +m' 2 )/2. ma ^g . is a polynomial function 
of the u(2) linear invariant N\ + N2 given by: 

(2i + l)! 



{mim2} g J (iV 1 + Ar 2 ) = 



\Jmax H~ J 71 



i)! 



-j] 



(13) 



where A^ = X x (X - 1) x ... x (X - k + 1). 

We note that when j — j max , the operator IT!?]) reduces to 

the identity and 



[mi -m 2 ]q'(jmax) — [mi -m 2 ] yi ( jmax ) 



(14) 



the minimal covariant component of which is given by (J6j) 
All phase conventions have been settled so that under her 
mitian conjugation (f ) and time reversal (K-t) we have 

[mi -m 2 ]^"t(j) = (—ly-m ^[mi -m 2 ]^"0) ^ 

_ ,-mi— m 2 [m 2 -mi]<r(j) 
1 - 1 m I 



(15) 



2.1.2. Matrix elements within the standard basis 

With the Wigner-Eckart's theorem we have in terms of 
su{2) Clebsch-Gordan coefficients (CG) [30j: 

(([n"0}j"m"\^ - m ^\[n'0\j'm')) = (2f + I)' 1 / 2 

xC m m ^* ([n''0}f\\^- m ^T^\\[n'0]f) 
(J f) m" V ' 

= p m! ([n"0]j")* 

([m[- m ' 2 }j [n'0]?) 

x(2j" + I)" 1 ' 2 ([n"0}j"\\^- m ^T^\\[n'0}f) , (17) 



where * denotes complex conjugation. The notation for 
the F symbols, which retains the full u(2) 2> su(2) labels, 
is useful when symmetry adaptation in a point group is 
performed. Reduced matrix elements (rme) for all opera- 
tors are obtained with 



(18) 



[n"0]i"|| [mi - m2] r«||[n'0]/J = «5„» n / + mi- m2 
(2j + l)(n' + m' 1 + l)\(n' - m' 2 )V 



jq [mi - m 2]q-{3) j(^ t ~ 1 — f— ^mi[mi -maifU) 4 



(16) 



(n' - m 2 )! [ (jmax + j + 1)! (jmax ~ j)\ 

2.2. Tensors adapted to an u(3) dynamical algebra 

For the applications we have in mind, the initial as- 
sumptions are the following: 

• A molecule, with point group symmetry G (G = C2 V ), 
admits in its full vibrational representation two non de- 
generate modes with close enough frequencies. 

• The interaction of these modes with other vibrational 
modes is sufficiently low so that in first approximation it 
can be neglected. As a consequence a separate study tak- 
ing into account the degrees of freedom associated with 
these modes only is possible. 

• This study, made within the frame of an u(p+ 1) dynam- 
ical approach, requires the introduction of a dynamical or 
non-invariance u(3) algebra to which we associate the ele- 
mentary boson operators {bf , &i}j=i,2,3- 

2.2.1. The algebraic chain u(3) D u{2) D su(2) D so{2) 

^From the preceding assumptions the space of states is 
a carrier space for the so-called totally symmetric (or most 



degenerate) irrep [NQQ] — [N 6] o f u( 3) which subduces to arbitrary covariant operator is obtained through: 



[nO] (n = 0, 1, • • ■ TV) in w(2) [29|, [3l(. Then all operators 
which may act within this irrep are of symmetry [z, 0, — z] 
with z = 0, 1, . . . , N (see for more details). 

Those which are maximal in w(3) have the form: 

(iV-z)P 1/2 



T 



TV! z\ 



(19) 



z - z) 
(max c ) 

where the a 2j jv coefficient is determined through the fol- 



lowing normalization condition: 



< n,0,N -n\ T 



[zO-z] 
(max c ) 



\Q,Q,N>=8 Z 



(20) 



The notation \rii, ri2, n 3 ) for the states is that of the u(3) 
canonical chain. |0, 0, iV) and \n, 0,N — n) represent re- 
spectively the state with zero excitation quantum and the 
state with n excitation quanta maximal in u(2). Then 
the semi-maximal operators of the u(3) dynamical algebra 
write as: 



T 



z -z 
mi ~m 2 
rrii 



= G{z,m 1 ,m 2 ) b+ mi b 2 m \ 



(21) 



where G(z,mi,m 2 ) is an operator valued function invari- 
ant in it (2) and defined by: 



g(z,m 1 ,m 2 ) 



mij \m 2 
(z + m 2 + 1)! (z + mi + 1)! (N - z)\ 



(mi+m 2 + l)!(2z + l)!iV!z! 

It (-l) t+m * ( z ~ mi \ (rai+m 2 + l)! 
lt=0 V t J {mi+m 2 + l + i 

v {Z - m2)l A N 1+ N 2 -m^ 



t)\ 



(22) 



(z — TO2 - t)\ 

x (N 3 - z + mi + u)^} b+ z ~ mi - u b 3 

with u — inf(z — m\,z — m 2 ). From equations (0 1141) 
it appears the left member of equation (I2ip is, within a 
phase factor, the minimal covariant component of an ITO 
within the su(2) D so(2) chain. So from the results in 
section f2.1.1( and taking into account that C/(Z, mi,m 2 ) 
commutes with the su(2) ladder operators J+ and J_ , an 



T 



z -z 
mi -m 2 ] (j) = (-l) m2 [mi!m 2 
m 

x g(z,m 1 ,m 2 ) 1™-™>1t<&, 



(23) 



where the various labels may take the values < z < N, 
< mi < z, < m 2 < z, j ' = (mi + m 2 )/2, — j < m < j. 
The phase of the preceding operators have been chosen so 
that 



T 



z - z 

TOl - TO 2 ] (j) 

m 



t 



z -2 
m 2 - mi] (j) 



(24) 



/c t T 



z - z 
[mi -m 2 ] (j) 
m 

xT 



Kt 1 = (-i) mi 

z - z 
K -m 2 ] (j) 
m 



(25) 



2.2.2. Matrix elements 

The covariant states, associated with the representa- 
tion [iVO] of u(3) and adapted to the subduction u(3) D 
u(2) D su(2) D so(2), are denoted |[7V0][n0]jm)) in the 
following. They may be generated from the zero quantum 
excitation state | [AT 6] [0 0] 00)) = |0, 0, iV >, through the 
relation: 



\[Nb}[n0]jm)) = T 



n — n 
[n 0] (j) 
m 



\[N0][Q 0]00)). 



(26) 



The matrix elements for operators (|23p , computed in basis 
I, are then given by: 

z - z 

[mi - ma] (j) 

771 



(([iV 6] [n" Q]j"m"\T 



\[N0] [n'Ojfm')) 



i In m m ' (j")* 

(2j" + l)- 1 / 2 C U 7 

(j /) m" 



(27) 



[NO] [n"0]j"\\T 



\\[N0] [n'0]f 



z - z 
[mi - m 2 ] (j) 

Using the expanded form (|22p of Q(z,mi,m2) we obtain 
the rme in the form 

z -2 
[mi - m 2 ] (j) 



([iVO] [n"Q]j"||T 



1 1 [TV 6] [n'O]/ 



z \ / z \ (z + m 2 + 1)! (z + mi + 1)1 (N - z)\ 



miJ\m 2 J (mi + m 2 + 1)! (2z + 1)!7V! z! 

mi! m 2 ! (iV - n')! (N - n' + m 2 - mi)! 1 1/2 
(JV - n' - z + m 2 + u)\ (N - ri - z + m 2 + u)\_ 

t (z — niA (mi + m 2 + l)!(z - m 2 )! 



(mi+m 2 + l + t)! 

(n' — m 2 )!(iV — n' + m 2 — z + u)\ 
(z — m 2 — t)!(n' — m 2 — t)\(N — n' + m 2 — z + t)\ 

x([N0] [n"0]f \\l m ^- m ^T^\\ [n'0]j'), (28) 

where the remaining rme is given by (|18p with m'i = mi, 
m' 2 = m 2 . 



3. Application to XY2 molecules 

We make the assumption that an appropriate dynami- 
cal algebra for the description of the vibrational spectrum 
of these molecules is u(3)s x u(2)b, where the indices S 
and B respectively refer to stretching and bending modes. 
Below, we first specify the notations which are used in for 
each type of mode and next perform the symmetry adap- 
tation in C 2v for states and operators. 



3.1. Bending mode v 2 

3.1.1. Standard states and operators 

For the Schwinger's realization of su{2)b we take ©: 

^J+ = b+b 5 , ( fc )j_ = &+6 4 

W J z = ±(7V 4 - N 5 ) = ±(b+b 4 - b+b 5 ), 

and (b) l[ 2) = N 4 + N 5 = N b for the u(2) B linear invariant. 
The standard covariant basis © is written (N b = 2J b ): 



\[N b 0}J b m» = (-l) 2 ^+™ \J b -m,J b + m> 



r_i\2J b +m 



[(J b + m)\(J b -m)\] V2 4 
We may also set 

J b — m = n 4 = v 2 
J b + m = n 5 = N b — v 2 



L+J b ~mi+J b +m 



0,0 > 



(29) 



Jb 



n 4 + n 5 N b 



n 5 - n 4 _ Nh 
2 ~~ 2 



v 2 . 



Thus all results of section [2~T1 can be used with the appro- 
priate change of indices. However as N b is associated with 
the maximal number of bending states and that only oper- 
ators acting within the [N b 0] irrep of u{2) are allowed we 
necessarily have m 4 = for the general operators built 
before. The vibrational operators are thus (0 < j b < m 4 ) 



[m 4 ~m 4 ]g(j b ) _ p b -m 4 



[m 4 0] T (^) x [0-m 4 ] T (^) 



= {m4 ' m4] g Jb {N 4 + N 5 ) K-"*«]bW*), (30) 
with m.4 = m' 5 = j b and 



{m 4: m 4 }g (7V 4 + iV5) 



(2J6 + 1)! 



(m 4 +j b + 1)! (m 4 - j b )\ 
x (jV 4 + Ar 5 -i 6 )[ m ^. 



1/2 



(31) 



Their expanded expression in normal ordered form is: 



[m 4 -m 4 ]g(j b ) _ p b —m 4 ^2qi 



V 2 



V 2 



/m 4 \ /m 4 . 



' (^ ^) m 



b+^- qi bi^ +qi b A ^ +q2 & fi 



(32) 



5 



With equations (fT71 [T8| their matrix elements are 

({[N b 0}J b m / '\^- m ^B^\[N b 0}J b m')) = (2J b + l)~l 
m m' (\N b 0]J b )* 

xF 

([m' 4 -m' 4 ]j b [N b 0]J b ) m" 
x {[N b 0}J b \\^- m ^B^\\[N b 0}J b ^ , (33) 

with for the rme 



([N b Q]J b \\^- m ^B^\\[N b Q]J b ) = 



-3b 



{N b -m 4 )\ 
(2j b + l)(N b +j b + l)\(N b -j b )\- ' 



(34) 



(m 4 + j b + 1)! (m 4 - j 6 )! 

3.1.2. Symmetrization in Ci v 

For the mode with A\ symmetry there is a priori 
no difficulty. All operators I"^ -m^gCA) are of A\ species 
as well as any linear combination 



[1714 —7714] 



Ktib) \ , [m4-m4]rim [m 4 — m 4 ] o(jb) 



(35) 



where I 1 ™ 4 _r " 4 lG is a unitary matrix and r = 1, • • ■ , 2j b +l 
a multiplicity index. However a consistent tensor formal- 
ism imposes some restrictions upon the G matrix. 
The matrices for the irrep of Ci v , which reduce to char- 
acters, being real, the metric tensor may be chosen (and 
is usually chosen) identical to the identity. This traduces 
by: 



for any symmetrized operator, and in particular for the 
elementary ones 

6 (AOt = \(M) _ 

For our problem we will not undertake a general discus- 
sion as was made for E modes .241] . For i>2 we only need 
the similarity transformations \- Nb °>G for the states and 

[m 4 -m,i]Q f or the operators. Also from the expression (j3"Tj) 
f or {m 4 ,m 4 }g (JVlt+iVs), with Ai symmetry and invariant 
upon time reversal, we see that we only need to determine 

[m 4 -m 4 ] G for the operators I™ 4 -^SC^) (j 6 = TO4 ). That 



is (equation (|30l) ): 

[m 4 -m 4 ]g(jb) _ {m 4 ,m 4 }g (jV^-fJVs) 

x y ' [m' 4 —m' 4 ] Qm^ [m' 4 -m' 4 ] g 



(jb) 



(36) 



• Symmetrized states. From equations (|9I16[) we have 
AC t |[JV 6 0]J 6 m)) -ACt^ '^)^- 1 ^ |0,0> 



l)^\[N b 0]J b m)). 



(37) 



If we set 



|[JV 6 0]J 6 rAi)) = e ie |[A^0]J 6 m)), 

we get symmetrized states invariant upon time reversal if 
e 2ie _ ^ — 1)^6 = i 2Ar *>, thus a possible choice is: 



\[N b Q]J b rA x )) = \[N h Q\J h mAi)) = i" b \[N b 0]J b m)) 
= {-\) J »- m |n 4 = J b - 171,115 = J b + m>, 

which amounts to choose 



[m 4 0]s~<m _ ■m 4 sr 



(38) 



(39) 



• Symmetrized operators. Likewise from the properties 
(|15l I16|) under hermitian conjugation and time reversal of 
the standard tensors, it may be shown that we may built 
the hermitian symmetrized operators 

[m 4 — m 4 ] ofe) 2. f[ m i — ™u] PZ(n) 

d |to|eAi ~~ ^2 V ™ 

+ (— 1) £ (— iy b+m I™ 4 -^IgO't)^ (40) 

[m 4 — m 4 ] KJU&J -rn 4 [m 4 — m 4 ] fc}(ii>) 

D 0Ai ~ 1 a ' 

with e = 0,1. Operators characterized by e = are in- 
variant upon time reversal (resp. non invariant upon time 
reversal) for j b even (resp. j b odd); it is the reverse for 
those characterized by s = 1. Operators "diagonal in D2" 
[m 4 -BMlgW are i nvar j an t upon time reversal. We may 
thus write an effective bending Hamiltonian: 



Nb m 4 
m 4 =0j' b =0 



(41) 



• Matrix elements in the symmetrized basis. Several meth- 
ods can be used for their computation. The simplest is to 



use the general formalism, that is to perform the change of 
basis associated with the symmetrization process of states 
and operators in equation (|33p . For the general operators 
we obtain 



with 



«[«iO]J t m"A 1 |[""-""lBg||[«iO]J l mM,)) = 

{2Ji+lr , F rM m'A, (Wl* 

([m'^-m^jb [N b 0}J b ) m"Ai 

x {[N b 0}J b \\^- m ^B^\\[N b 0}J b ^ , (42) 

where the rme are those of equation (|34p ; r stands for 
r± = |rn|e or according to the case and the symmetry 
adapted CG coefficients are given by: 

rA x m'Ax ([N b 0]J b ) 



F 



([m' A - m'Aj b [N b 0]J b ) m"A x 



m.m' ,m' 



[m 4 — m 4 ] /im * [A'';, 0] /-tm * [N^ 0] 



xF 



E 



TO to' ([JV 6 0]J t ) 

([mi-mi]j6 [iV 6 0]J 6 ) to" 



xF 



m 



m 



(43) 



([JVfcO]J 6 ) 
{[m' 4 -m' 4 }j b [N b 0}J b ) to" 

where the last equality follows from (f39|) . The analytical 
expressions for these coefficients together with those of the 
various matrix elements are given in Appendix \Xi 

3. 2. Stretching modes v\ , V3 

3.2.1. Standard states and operators 

For the su{2)g Schwinger's realization we take: 

= b+b 2 , (s) J- = b+h 
{s) J z = \ (N 1 -N 2 ) = ^ (6+6! - 
where the indices i = 1, 2 are linked to the two bonds. 
With regard to the notations of section 12.21 we make the 
substitutions N — > N s , j — * j s for the states and T — > S 
for the operators. We have then with (j2"6f (j s = n s /2): 

n s — n s 
[». 0] (j a ) 



^2) 



(44) 



|[JV,6][n.0]j.m)> = S 



|[JV.0][00]00», (45) 



|[iV s 6][0 0]00» = |0,0,JV a ) - (NJ)-V 2 b+ N °\0,0,0). 

We may also express these states in various forms replacing 
the S operator in 



\[N s b][n s 0]j s m)) = (n s \)^ 2 g(n s ,n s ,0) ^S^}\0,0)\N s ) 
= (n s !) 1/2 G(n s ,n Sl 0) \[n s 0]j s m))\N s ) 
= (-1)*+"^ \ j.-m,j. + m,N„-n B ). (46) 

We thus have the correspondence: 

n s n s 

— - to = ?ii , — +ni = n 2 , N s -n 8 = n 3 . (47) 

Arbitrary vibrational operators in the algebraic standard 
u(3) D u(2) D su(2) D so(2) chain are obtained from re- 
lations (l2"2l |2"3"|) and the results for the covariant operators 
[mi -to 2 ]j>u)_ ^11 their matrix elements are given by equa- 
tions (|2~T1 12"5|) . In particular, operators which are "diagonal 
in n s " are characterized by mi = m-i. 

3.2.2. Symmetrization in C 2v 

The computations have been made with the conven- 
tions given in Table [TJ The indices i = 1, 2 being associ- 

Table 1: C 2 v character tabic 



Civ 


E 


C 2 {Oz) 


a v (xz) 


a v (yz) 




1 


1 


1 


1 


A 2 


1 


1 


-1 


-1 


Si 


1 


-1 


1 


-1 


B 2 


1 


-1 


-1 


1 



ated with the two bonds we have: 

Pc 2 h Pel =b 2 ; Pa h P- 1 = b 2 ; P a , h P-, 1 = b h (48) 

where we set a = a v (yz) and a' = a v (xz). Hence the set 
(bi, b 2 ) (or (6 + , b 2 )) span a representation A\ +B\ of C 2v . 



• Symmetrized operators. The relations in (|48p allow first 
to determine the transformation laws of the standard ten- 
sors (i = (mi +m 2 )/2): 



[mi -rn 2 ]c<{j) p—1 _ ^_^y m i [ m i — maJjj'O') 



(49) 



where i? = C 2 , cr; they are obviously invariant with respect 
to cr'. We may thus build the symmetry adapted tensors 
(m > 0) 

[mi -ma] _ 0(mi, m 2 , |m|, T) [mi _„ l2] y) 

V2 (50) 

+(-l) mi+e[mi - m2l S°^), 



with e = for T = A x and e = 1 for T = B x . For m = 
we simply have 



[tm -wal^W) = g(mi , m2, 0, r) [rni ^ m2] , 



with 

T = Ai for mi and to 2 even 
r = B\ for m-i and to 2 odd. 



(51) 



(52) 



In equations (|50l l5Tj) 6* is a phase factor to be fixed next. 
In order to obtain a correct description (in terms of allowed 
symmetries in C 2 „) for the states associated with the irrep 
[NO] of u(3), we must impose that (or 63) belongs to 
the A\ scalar representation of C 2tI . With these conven- 
tions, it appears that the Q(z, mi, m 2 ) term , given by (|22[) . 
is invariant in C 2l) . Consequently the general operators 
(|23|) transform, under the action of the C 2l) generators, as 
the standard operators (relation (|4T)|) ). This property, al- 
lows to symmetrize both type of operators with the same 
orientation matrix and we set for the u(3) symmetrized 
operators: 



s 



z -z 
[mi - m 2 ] (j) 

Imir 



E[mi — m 2 ] /nr m Q 
^Imir 



z - z 
[mi - ma] (j) 
m 



(53) 



where the sum is limited to the values m = ±|m| (to 7^ 0). 
The choices for the phase factors have been made so that 
the operators ([53"|) satisfy: 
t 



S 



z - z 
[mi - to 2 ] (j) 
\m\T 



S 



z - z 
[m 2 - mi] (j) 

Imir 



(54) 



Under time reversal they are unchanged when mi or m 2 
are zero and multiplied by (— l)-?+l m l+ E when mi,m 2 =/= 0. 
This leads to 



[mi "" l2l G|™|r - "4 0(mi>m 2 , |m|, T) , m > 



V2 



[mi -m 2 ] G ^m = e(mi)TO2) |m|, T)(-l) 



mi+e 



V2 



(55) 



[mi ~ m2lG lo°ir = ( m i' TO 2,O,r) ,m = 



with 



0(mi,m 2 , |to|,T) = i mi f+\ m \+ e toi,to 2 ^ 
6»(toi,0, |to|,T) =i mi 

0(o,TO 2 ,H,r) = i m - (-i)i" l i+ e . 



(56) 



With the results in this section, keeping terms which are 
diagonal in n s only, we may write the effective stretching 
Hamiltonian: 



z=0 mi=0 |to s | 



z - z 
mi - m x ] (j s ) 
\m s \ Ai 



(57) 



with j s = toi and where the sum S* is over |m s | values 

such that j s + \m s \ = mi + \m s \ be even. 

• Symmetrized states. These are obtained with equations 

(251 1351 B71 153): 



|[JV,6][n.0]j.|m|r)) = S 



n s — n s 
[n s 0] (js) 

imir 



|[AT S 6][0 0]00», 



(58) 



with [ n >°]G obtained from flHHl OS}. Explicitly we have for 
the various types of local states: 
Local states {ni,n 2 } n\ = n 2 . 



They are associated with m = (Eq. l47l) which implies j 
integer or (as j s = n s /2) n s even. 

, n s n 



|[iV s 0][n fl 0b- fl 0Ai)) = (-1)- 



,N s -n s ). 



(59) 



Local states {ni,ri2} n\ ^ n 2 . 

We have then \m\ ^ and still setting e = for V = A\ 
and £ = 1 for T = Bx: 



|[jV s 6][n s 0]j s |m|r)) = 
s/2 



— {It -m, — + m,N s 



-(-iy\f + m,f-m,N s -n s )}. 



(60) 



With our phase convention they are all invariant under 
time reversal. 

• Matrix elements in the symmetrized basis. They are ob- 
tained with a method similar to that used for the bending 
mode. From equation (|27| the transformation to sym- 
metrized states and operators gives 
z -z 
[mi - ma] (j) I*')) = <->/' + D ' ' 

Imir 



«*"|S 



xF 



[iV s 6]K'o]^||S 



\\[N s 6][n' s 0)f s (61) 



|m|r \m'\T' (K'0]j«)* 

([mx-malj [<0]ii) |m"|r" 

z - z 
[mi -m 2 ] (j) 

where we set |\&)) = \ [N S 0] [n s 0]j s |m|r)) and the rme are 
given by (|28|l with the appropriate label substitutions. The 
symmetry adapted CG coefficients are obtained with 

p \m\r Kir" (K'oK') = 
(K-m 2 ]j kok) | m "|r" 

w \ ^ [mi —1712] s~im * \n' 01 f~ifn * [n'' 0] rfm" 

x 2^ u |m|r U |m'|r' ^im'^r" 



m.m' ,m" 



x F 



(62) 



m' (K'O]^') 
([mi-m 2 ]j [n' s 0]^) m" 
with j = (mi + m 2 )/2, ^ = n' s /2 and j^' = n"/2. Also it 
is important to note that these coefficients are a priori de- 
fined only for n" s = n^+mi — m 2 . The matrices I™ 1 _m2 ]G 
and '" s °1 G are given by equations ([531 IS"6")) . The analytical 
expressions for these coefficients are given in Appendix [B] 



3.3. Stretch-bend interactions 

Results in sections 13.11 et 13.21 determine all operators 
adapted to the study of isolated bending and stretching 
modes and which may appear in the Hamiltonian or tran- 
sition moments. 

Taking into account stretch-bend interactions introduces 
coupling operators which may be formally written 



(?s b ) = [5 cr.) x B (r fc)] 



(Fa!,) 



For our problem = A\ and the CG for the C21, group 
are trivial ; we thus simply have 



O Cr.) =5 <r.) B (A0 



(63) 



with T s = A\ for Hamiltonian terms. The operators ([63| 
may be written in various manners depending on the case 
and also depending on what we mean to represent, for 
instance operators in the untransformed Hamiltonian or 
effective ones. 

3.3.1. Method in the case of a polyad structure 

The chosen dynamical algebra assumes u>\ ~ uj^(= lu s ) 
and we have to take into account in the effective Hamil- 
tonian the resonance with the bending mode which deter- 
mines the polyad structure. We assume 

uj\ ~ LJ3 = lu s ~ k ui2 = kcub- 

Within our formalism the operator N% + N% = h s , with 
eigenvalue n s , represents the "number of quanta of stretch- 
ing" operator ; N4 with eigenvalue 714 = rij the "number 
of quanta of bending" operator with 
N 4 + N 5 N 4 - N 5 N b 



N 4 = h b 



2 2 + Jz 



(64) 



To a given P polyad we may associate the P operator 

P = n b + kn s = N i + k{N l + N 2 ), (65) 

which may be expressed in terms of the ITO defined pre- 
viously knowing that: 

N 4 = V2(^4°l+^Bil); 



N!+N 2 = -{N s 



3N, 



S 



10 - 1 
[0 0] (0) 



}■ 



An O operator which conserves the P quantum number 
associated with a given polyad must satisfy the condition 
[K, O] = 0, be of species Ai and invariant upon time re- 
versal if it belongs to the Hamiltonian expansion. 
To determine the possible O operators it appears that it 
is better to work first in the standard algebraic chain 



«(3)s 
[N s 6] 



u(2) B 
[N b 0] 



u(2) s 
[n s 0] 



u(2) B D 
[N b 0] 



su{2) s ® su{2) B D so{2) s ® so(2) B 
j s = n s /2 j b = N b /2 m s m b 

where the indicated symmetries are those of the states. 
We thus start from the operator basis 



SB = S 



z -z 
[mi - m 2 ] (j s ) 
m. 



[m 4 —mi] K?(it) 



(66) 



with j b = m' 4 and j s = (mi + m,2)/2. A straightforward 
calculation gives 

[P,SB] = [-rat + fc(mi - m 2 )]SB, 

thus the condition [P, SB] = is satisfied if 



m b = k (mi — m 2 ), 



(67) 



with k = 2 for the considered XY% molecules. This con- 
dition being independent of m s , it appears that, in order 
to determine hermitian interaction operators having also 
a determined behavior upon time reversal, it is preferable 
to keep the standard form for the bending operators and 
to take symmetrized operators for the stretching ones. We 
mainly have two cases: 
• m b — then mi = m 2 

In this case we can also take directly symmetrized bending 



operators, which gives the hermitian operators 



o sb (\) = S 



z — z 
[mi - mi] (j s ) 
|m s |r s 



[mi —mi] n{jb) 



(68) 



satisfying upon time reversal (see sections (13.1.21 
JCOs^JC- 1 = (-l) mi +l^l+ e O sb (l). 
• m b 7^ 

^,From the set (|66j) we may build the hermitian operators 
(equations fTSl lM jl ) 



O sb (2) = -j= < S 



z -z 
[mi - ma] (js) 



[m 4 -mi]n{j b ) 



+ (-i) T (-i)A+ m » s 



[m4-t7u] B CA) I ; ( 69 ) 



z - 2 
[m 2 - mi] (j s ) 

|m s |r s 

with r = 0, 1 and where mj > and mi > mi is assumed. 
With the properties established in sections (|3.1.2[ I3.2.2p ) 
one shows that 



K.O sb (2)K- y = 



(-l) r (-l)A(-iy'+l m 'l +B O a6 (2) (a) 



(-iy(-iy>o sb (2) (b) 

where case (a) (resp. (b)) is for m 2 7^ (resp. m 2 = 0). 

3.3.2. First operators in the Hamiltonian for k = 2 

Since the operators B only differ from the B ones through 
a constant function within the irrep [N b 0], we can make 
in equations (f6"51 16"5)) the substitutions 

[mi-mi]n(j b ) [m' A — m' 4 ] p>(j b ) /■ _ i\ 

to define the terms (and the parameters) of the effective 
Hamiltonian. 

• For m b = (and 7^ 0, z 7^ 0) the operators (|68[) . with 
r s = Ai and (— l) m i+l m »l+ E — 1 are products of operators 
belonging to Hb and Hs- 

• The cases m b 7^ correspond to non trivial (that is 
non diagonal) interaction operators. With k = 2 and as 
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mi and m 2 are integers the even values of mj, alone are 
allowed in (|69[) . with a minimum value rribmin — 2. As a 
result m' imm = j hmin = 2 and m imi n = 2. We have then 

(mi -m 2 ) 

min — 1 from which we deduce j s m in — 

1/2 and 

|m sm i„| = 1/2. We thus have the first allowed values for 

the labels in the operators (G9): 

i) mi = 1, m 2 = 0, j a = 1/2, \m s \ = |1/2| ^ r = 

zzj mi = 2, m 2 = 1, j s = 3/2, |m s | = \l/2\ => r = 

m) mi = 2, m 2 = 1, j s = 3/2, |m s | = |3/2| =>• r = 1 

We thus obtain: 



0.i(2)i 



1 

71 





z - z 




[1 o] (|) 







(70) 



with z = 1, 2 ■ • ■ . The operators with z > 1 are anhar- 
monicity corrections to those obtained for z = 1. 
We find next: 



O sb (2) 2 




(71) 



[2-2] B (2) 



h.c. 



(72) 



z - z 

a b (2) 3 = 4= <! S [2 - 1] (i) 

h 

with z = 2, 3 • • • . 

We also have other possible operators [^"H^W with 
m-4 > 2 in (|70M72p ; we restrict here to those which are of 
lowest degrees in creation and annihilation operators and 
which will be used in the next section. 

4. Application to D?S 

4-.1. Effective stretching Hamiltonian 

Up to the second order, that is for z < 2, the general 
expansion (|57|) involves terms with the following values for 



the labels 

z = mi = j s = 

z = 1 mi =0 j s = 

mi = 1 j s = 1 

z = 2 mi =0 js = 

mi = 1 j s = 1 

mi = 2 js = 2 

mi = 2 js = 2 



|m s 
|m s 
|m s 
|m s 
|m s 
I TO, 
|m s 



= 
= 
= 1 
= 

= 1 

= 
= 2 



With the results of sections [2 . 2 . 1 [ |3~2~21 Hs can be written 
in terms of elementary boson operators as 



^-[Nz-\{N 1 + N 2 )] 



+ < 2 >4 0) ^gf {3A 3 (A 3 - 1) - 6(iVi + N 2 )N 3 
v 30 



+ (Ai+A 2 )(Ai+A 2 -l)} 

+ {2} iT } ^# {4A 3 - (iVi + N 2 - l)}(b+b 2 + b+h) 
V 10 



(73) 



+ ^ (6+ 2 fe 2 + b\%\ ib+bihh) 

+ ^^(bfbl+bt%l), 



with a 2 ,Ar s = [2Ns(N s - l)] - *. 

For the fitting procedure, it can be rewritten to the more 
convenient form: 

H s = a (N! + N 2 ) + ai(A 2 + A 2 ) + a 2 AiA 2 + a 3 Y^ 

+a i {N 1 + N 2 )Y {Al) + a 5 [Y {Al) x Y {Al) ] {Al) , (74) 

with Y^ Al ^ = b^b 2 + b^bi and where we removed the part 
of Hs which depends upon the operator N s only since 
the latter takes a constant value within the irrep [A 0] of 
u(3)s- Also we set: 

f 2 }^(2) «2,jV _ { 2 } ? (2) «2,jV 

V6 2 V2 ' 



a 2,iV s 



1 - a, w \Wi l0 \ — + < 2 >t (2) — 1 
1 — «2,iVsL V 3 1/5 



a 2 = « 2 ,^y|[{ 2 >4 0) 2V5 - ^2 ~ {2} t 2 2) V3], 



a 3 



Qf4 



(i)f(D. 



r _ + ™ f <V K5 -j(4 W . + l,, 



_{2}r(i) 



*i '(X2,N S \/ ^ 



Q'5 



= {2} ? 2 2) «2,.V S — 
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4-2. Effective bending Hamiltonian 

The general expansion (|41|) may also be written, using 
the results in section (13.1. 2ft : 

m 4 =0j b =0 

with jb = m' A and {™ l 4,«u}g^ defined in (|3"Tj) . As a given 
jb (or m! 4 ) value appears for all 777,4 > jb values and since 
N 4 + N5 takes a constant Nb value within the irrep [Nb 0] 
of u(2)b we may set 

n B = J2 ^ K ~ <]b oaI > ( 75 ) 
with the effective parameters: 



^0&) _ {m,4}^(if>) {m 4 ,m 4 } f 

For instance up to second order we have: 



■1/(2) _ r(0) r(i) [1 -1] R (l) , r(2) [2 -2] R (2) 



0A\ ^ % 



(76) 



(77) 



7(1) 



+ l JL=(N 4 ~N 5 ) 



V2 



f(2) 

' I, 



+^[bt 2 bi+b+ 2 bi-4b+b+b 4 b 5 }, 



(78) 



where the last form has a clearer physical meaning, with 
parameters (3$, /3 4 given by: 

^.5. Effective stretch-bend interaction operators 

As noted before these interaction operators may be di- 
vided in two groups. In the first one we have products of 
stretching and bending operators diagonal with respect to 
n s and rib- In the second one are those obtained from the 
properties of the w(3)s and u{2)b dynamical algebras and 
which take into account the approximate resonance be- 
tween the stretching an bending modes. Keeping only op- 
erators of lowest degree the stretch-bend effective Hamil- 
tonian can be expressed as (equations l[6"8]). ([70))): 



= 72 (b+b 3 bfbl + b+b 3 b+ 2 bl + h.c.) = 72 (80) 

It is worth to analyze more deeply this last operator which 
indeed traduces the resonance vi{Ai) ~ v^{Bi) ~ 2v-i{A\). 
Thus the Hamiltonian matrix, already divided into two A\ 
and B\ blocks , is subdivided into sub-blocks characterized 
by the polyad number P — 2(n\ + 712) + nb. Within each 
P-block, Hgg connects states which are not diagonal nei- 
ther in n s nor in nb- 

O 1 ^ can also be written (bf + b\ )bl[bf 2 b 3 ] + h.c. In 
this form it is clear that the dependence of its matrix ele- 
ments upon the quantum numbers m, ni and nb are simi- 
lar to that of a usual Fermi interaction operator. However 
contrarily to the later, which leads to convergence prob- 
lems for high values of the quantum numbers, the other 
factor the matrix elements of which behave roughly as 
(N s — n s ) x / 2 (Nb — rib), has a damping effect as n s and 
(or) rib increase (N s and Nb fixed). 

We already defined a similar operator adapted to XY3 
pyramidal molecules in 



34 



351 ] and already proved that 



this operator does not require the knowledge of N s and 
Nb for low values of the quantum numbers n\, 77,2 and 77,3 
as it physically must be near the minimum of the potential 
function. 

4-4- Numerical application 

To illustrate the efficiency of our formalism, we apply 
it to the deuterate hydrogen sulfide. The hydrogen sul- 
fide molecule and its isotopic species are of interest for 
terrestrial atmospheric pollutant measurements. They are 
involved in the study of planet atmospheres and appear in 
the analysis of interstellar medium. Also, hydrogen sulfide 
is a good candidate to apply local mode models. We re- 
strict here to the D2S molecule but will present a compar- 
ative analysis with other XY2 molecules in a next paper. 
Renaming, for simplicity, the parameters a$ (i=0, ...,7), 



H'sb = lo(Ni + N 2 )N 4 + 7l 7V 4 



(79) 
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the best fitted Hamiltonian is defined as follows: 



H = ao(2Vi + N a ) + ai(JV? + Aff) + a 2 iViiV 2 + a z Y^ 



+ a 4 N 4 + a 5 N( + a 6 Oi/ B + a r (Ni + N 2 )N, 



(81) 



However, taking into account the N s and N b quantum 
numbers, there are 10 parameters to determine in our 
model since the value of these numbers appears explicitly 
in the matrix elements of To obtain the optimum 

values of these quantum numbers, we simply operate as fol- 
lows. An initial set for the parameters 4°^ (i = 0,1,4,5) 
is obtained from the experimental values of the lowest vi- 
brational bands. This allows to have approximate values 
for the stretching levels when n s quanta are localized on 
one bond and for the bending levels: 

£j^'(n s )—a n s + a 1 n s , & b (n b j — a 4 n b + a 5 n b . 

As for any anharmonic potential (Morse or modified Posch- 
Teller for instance) the n max value of the vibrational quan- 
tum number is given by the extremum of the E^ (n s ) and 
E b (n b ) curves. With 



i (0) 
,(°) 



1927.5855 cm" 



4 0) 



-24.4279 cm" 



4 J = 858.2604 cm- 1 , 4 = -2.8564 cm" 1 , 



we obtained in this way: 



150.23, 



«o 0) , - 4 0) 

T^s max ' 77v\" 39.45 , ^bmax TTyT 

2af> 24 ' 

hence we could take (iVj 0) , N^) = (39, 150). However this 
high value for N b is not reasonable for the following rea- 
sons. As it is commonly accepted N s may be interpreted 
as the number of excitation quanta which, when concen- 
trated on one bond, may dissociate the molecule. Within 
a polyad P to which a |7V S , 0, 0} stretching state belongs, 
we also have states |0, 0, N S )\2N S , N b — 2N S ) which implies 
that N b > 2N Sl but on the other hand if we took A^ (0) 
much greater than 2ivi°^ , we would have many dissoci- 
ating states within polyads with P > 2ivi°' ) . It would 
be quite unrealistic to pretend that our model is capable 
to reproduce isolated stable states within the continuum. 



Another method to obtain a reasonable value for is 
to use for E^°'(n s ) the known experimental value of the 
dissociation energy Q: 

4% s + a { ° ] n 2 s ~ 32050 ± 50 cm" 1 , 

which leads to TV] -* = 24. Various fits were performed 
starting thus with (A s (0) , iV b (0) ) = (39,78) while the other 
a, (i = 0, • • • ,7) parameters were determined through a 
usual non-linear least square fit method. We noticed that, 
except for the a@ parameter associated with the ma- 
trix elements of which depend strongly upon the N s and 
N b values, other parameters remained almost unchanged 
(less than some percent of relative variation) while N s and 
Nb decreased. One of the indicator of the convergence of 
the fitting process was the minimization of the standard 
deviation 



<j{d,p) 



\ 



E 



(obs) 



where d and p are respectively the number of experimental 
data and the number of parameters included in the fit. 
It soon appeared that this indicator was rather insensi- 
tive to the N b value in a rather large range. A similar 
effect was already noticed in previous studies of pyrami- 
dal molecules 1351 ] . On the other hand using the initial 
(Ai 0) , iV b (0) ) = (24,48) values improved drastically the 
convergence to the minimum value c(22, 8) = 0.514 cm" 1 , 
thus reducing noticeably the computational time. The 
second indicator is, of course, the parameters stability at 
the end of the fitting procedure; the last variation Aa^ 
of the parameters fulfilled the condition |Aai/aj| < 10~ 7 , 
(0 < i < 7). The retained Hamiltonian (jSTj) led to the set 
of parameters given in Tabled] (these are given with values 
in parentheses which are la statistical confidence intervals 
in units of the last digits). 

Finally we compared the experimental dissociation energy 
with the value calculated from our model. This was done 
in two ways. First, removing all off-diagonal terms in 
the Hamiltonian (JSTJ, the energy of the |24,0, 0Ai) (or 
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Table 2: Parameters (in cm 1 ) fitted with 22 experimental data. 



a = 


1927.908(290) 


ai = 


-24.665(103) 


a-i = 


-0.845(288) 


£13 


-6.428(140) 


ff : — 


858.821(447) 


a 5 = 


-3.103(141) 


a& — 


0.005(001) 


a 7 = 


-10.488(150) 



1 24, 0, B{)) pure stretching state is computed with the an 
and ai parameters of Table [2] which leads to the dissocia- 
tion energy Ed — 32063 cm -1 , or more precisely if we take 
into account the parameter uncertainties 31996 cm -1 < 
Ed < 32129 cm -1 . The second way to calculate Ed is 
to diagonalize the Hamiltonian (|81| within the P — 48 
polyad with the parameters (i=0, ...,7) of Tabled The 
value obtained for the |24,0,0Ai), |24,0,0Bi) levels is 
E' D = 32062 cm -1 which is 1 cm -1 close to the previously 
calculated Ed value, which confirms the validity of our 
previous assumption to keep only diagonal operators to 
evaluate this dissociation energy. Taking again into ac- 
count the uncertainty on the parameters, we found that 
the diagonalization of the P = 48 Hamiltonian matrix 
with the two sets: 

{ai (min) } = {a^~ , a J" , , <Zg ,a§ , dg" , af } 
{a,i(max)} = {aj, af , aj, 03", 04 , % , a%~, ay"}, 

where a~ (resp. af) stands for the lowest (resp. highest) 
value of the ai parameter according to the la statistical 



confidence interval leads to 

E' Dmin = 31996 cm -1 < E' D < E' Dmax = 32128 cm" 1 . 

Both methods lead to similar values also consistent with 
the experimental results. 

Experimental and calculated energies of D2S vibrational 
levels are given in Table [3J The first column gives the 
normal notation {y\ v-i 1/3) of the level. The second one 
indicates the polyad P number. The third column indi- 
cates the levels in local notation with explicit symmet- 
ric (+) and antisymmetric (-) labels. The usual local 
mode notation has been adapted to our notation as fol- 
lows mn^, v « Hin^, = n&. Column 4 gives the eigen- 
values, whereas column 5 indicates the observed energy 
levels. Column 6 gives the difference Observed-Calculated 
energy. The last column shows that the eigenkets are close 
to the initial basis given in columns 1 and 3. We mention 
that all the experimental data used in the present paper 
are reported in [l^ . 

5. Conclusion 

We developed a formalism which allows a complete de- 
scription of vibrational modes in local XY% type molecules. 
Stretching, bending, interaction and transition operators 
have been built and analytical expressions for their matrix 
elements established in the chains u(2) D su(2) D so(2) 
and u(3) D it(2) D su{2) D so(2). Next a full sym- 
metry adaptation in the Civ molecular point group has 
been performed. This formalism has been applied to the 
D2S molecule where the 2:1 resonance between stretching 
and bending modes has been taken into account through 
an adapted Fermi-type operator the properties of which 
have been discussed. From a simplified model we derived 
reasonable values for the highest stretching N s and bend- 
ing Nb quantum numbers. Experimental data are repro- 
duced with a standard deviation a — 0.5 cm -1 and only 
8 effective spectroscopic parameters. The dissociation en- 
ergy calculated with these parameters is in good agree- 



mcnt with the experimental one and this also confirms our 
method for the determination of the N s value. Our ap- 
proach will be applied to other XY2 molecular systems in 
a next paper. 
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A. Symmetry adapted Clebsch-Gordan coefficients 
and matrix elements for 1/2 

The expressions for the ~ m ^G matrix elements are 
obtained from equations (I3U1 HU]) which lead to: 



m,0i 



^ImleAi — m >l m l 



(A.l) 



m, -m. 



As a result: 

• For r = we obtain the CG (03) 



F 



QAi m'Ax ([N b 0]J b ) 

(K - Kb'fe Wb0}J b ) m"Ax 

m' (J b ) 
= r m *C =i- m *5 m ,, m „, (A.2) 

(jb Jb) m" 
and for the matrix elements PS]) : 



(([N b 0]J b m"A 



II A I [m.4 — rm] isijb) 



x(2J b + l)--2 ([N b 0}J b \\^- m ^B^\\[N b 0}J b ) .(A.3) 



• For r = \m\s. With (|A.1|) we obtain in this case for the 

CG dm 



C 



ml A\ 


([N b 0]J b ) 


\ at, nl 7, '\ 

[IV b UJ Jb) 


Tli J±\ 


m ml 


(J») 


(jb Jb) 


m" 


-^ym , 4 -\-m £j 


—m ml 






( jb Jb) 



(A.4) 

This last expression is to be taken with m > (one can also 
set \m\ in the CG). We then have for the matrix elements 
62|): 



(([N b 0]J b m"A 



II A |[m 4 -m 4 ]n(jb) 



BCUWbOVbm'At)) 



1 i £ I rn ml (Jb) 
2J 6 + l)-'^= [C V b) 

V2 \ ( Jb J b ) m" 



+(-l) £ (-l) m * +m C 



-m ml (J b ) 



( jb Jb) ml' 
([N b 0}J b \\^~ m ^B^\\[N b 0}J b ^ , 



(A.5) 



where we took into account that the su(2) standard CG 
are real. 

B. Symmetry adapted Clebsch-Gordan coefficients 
for 1/1, ia, 

We give below the symmetry adapted CG coefficients 
used in the study of stretching modes. They are obtained 
from equations ([551 [561 16"2")) . The remaining standard su(2) 
coefficients are taken from [jjcj. 
• \m\ = \m'\ = \m"\ = 

This condition implies n' s — 2j' s , n" — 2j" and mi + rri2 
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or 



(K-m 2 ]j [n' a O]j' a ) OAi 



C 



o o (j'J) 



U f s ) o 

(n' s - m 2 )! (mi)! (m 2 )! 

(n' s +771! + 1)! 

(K + mi)/2)! 



= i 3 (n s + mi — + 1)- 



5 r,Ai- 



(B.l) 



(«-m 2 )/2)!(m 1 /2)!(m 2 /2)! 
As it is known the CG on the second line in non zero only 
if j + j' s + j'J i s even which corresponds to the selection 
rule r = A x . 

• |m| = |m'| = 0, |m"| ^ or |m| = \m"\ = 0, |to'| ^ or 
|m'| = \m"\ = 0, |m| 7^ 

It is easily checked that all F coefficients are zero in these 
cases. 

• |m| = 0, \m'\ ^ 0, \m"\ + 
or \m'\T> ([n'JO]j'j 



(K-m 2 ]j [ni 0]j' s ) \m"\T" 



x C 



m . \ m 



-|(B.2) 



(i is) V 2 

We note that the factor 1 + (— i) m i+ e + e just traduces 
the selection rule V x V = V" . In fact for the CG to 
be non zero we must have (— + e = 1, that is for 
mi even V — Ax and e' = e" that is V = T" (Ax or 
Bx) ; for m x odd, V = B x and e' ^ e" that is V ^ V" 

((r',r") = (A, #0 or (r',r") - 

• \m'\ = 0, |m| ^ 0, |m"| ^ 

In this case we necessarily have j' s integer so n' s even and 
r' = Ax. Here again a factor 1 + (— l) e+e appears in the 
computation, which is equivalent to T = T" . We find 

m\T OAx ([n'JO]j'J) 



F 



(K-m 2 ]j [n' a 0]j' a ) |m"|r" 

m (j^') 
x C <5|m|,|m"| <5r,r" 

(i is) m 



(B.3) 



• |m"| = 0, |m| ^ 0, |m'| 7^ 0. 

This case is similar to the preceding one. We necessarily 



have j'J integer so n" s even and r" = Ax- Here again a 
factor 1 + (— appears in the computation, which is 
equivalent to T — V . We obtain 



F 



|m|r 



\m'\T> ([n'JO]j'J) 



([ mi -m 2 ]j [n' s 0]f s ) 0A ± 

m — m (j'J) 
x C ^|m|,|m'| Sr.r 1 - 

U fs) 



jmi jj-\m\+e 



(B.4) 



We underline that in equations (|B.2i IB.3[ IB.4|) the su(2) 
CG in the right member is to be taken with m (or m') 
positive. 

• |m'| 7^ 0, |m'| 7^ 0, |m"| 7^ 0. 

In all cases a coefficient 1 + (— l) £+£ +£ appears in the 
calculation, which traduces the selection rule r" = T x T'. 
We obtain the following non-zero coefficients: 
- Tor \m"\ = Iml + Ito'I 



F 



\m\T 



W\T' ({n'jO] 3 'J) 



{{mx-m 2 ]j [n' s 0]j' s ) \m"\T" 

J_ m m' (j'J) 
V2 (j j' s ) rn + m' 

For m > m' > 



__ j-™2 i—j—\m\—e 



X^ C <5lm"|,|m| + |m'| <5r",rxT'- (B.5) 



|m|r 



Im'ir (K0] Js ) 



([mi-ma]j [<0]^) |m"|r'' 



1 m -m' (i s ') 
x — = C 



j' s ) m — to' 



F 



V2 (j 
For m' > to > 

|m|r |m'|r (K'o]j s ) 

([mi-ma]j K 0]j s ) \m"\T" 



<5|m"|,|m|-|m'| ^"Jxr'- (B.6) 



1 m —to' (j'J) 

X -7= C <5|m"|,|rn'|-|m| ^"Jxr'- (B.7) 

* U is) m — m' 



References 

[1] J. Lamouroux, L. Regalia-Jarlot, VI. G. Tyuterev, X. Thomas, 
P. Von der Heyden, S.A. Tashkun, Yu. Borkov, J. Mol. Spec- 
trosc. 250 (2008) 117-125. 

[2] Z. Zelinger, A. Perrin, M. Stfizi'k, P. Kubat, J. Mol. Spectrosc. 
249 (2008), 117-120. 

[3] T.A. Ford, J. Mol. Struct. 834-836 (2007) 30-41. 



16 



[4] A.J. Friedson, Icarus, 177 (2005) 1-17. 

[5] M. Smirnov, PhD thesis, IAO-Tomsk, Russie, LPM-Orsay, 
France 2001. 

[6] Y.C. Minh, W.M. Irvine, L. M. Ziurys, Astrophys. J. Lett. 345 
(1989) L63-L66. 

[7] E. Herbst, D.J. Defrees, W. Koch, Notices of the Roy. Astro- 
nomic. Soc. 237 (1989) 1057-1065. 

[8] R.J.A. Grim, J.M. Greenberg, Astron. Astrophys. 181 (1987) 
155-168. 

[9] J. Pliva, V. Spirko, D. Papousck, J. Mol. Spcctrosc. 23 (1975) 
331-342. 

[10] R.L. Cook, F.C. Dc Lucia, P. Helminger, J. Mol. Struct. 28 
(1967) 237-246. 

[11] J.R. Gillis, R.D. Blatherwick, F.S. Bonomo, J. Mol. Spcctrosc. 

114 (1985) 228-233. 
[12] C. Camy-Pcyrct, J.M. Flaud, A. N'Gom, Mol. Phys. 65 (1988) 

649-657. 

[13] I.N. Kozin, P. Jensen, J. Mol. Spcctrosc. 163 (1994) 483-509. 
[14] A.W. Liu, O.N. Ulcnikov, G.A. Onopcnko, O.V. Gromova, E.S. 

Bekhtereva. L. Wan, L.Y. Wan, S.M. Hu, J.M. Flaud, J. Mol. 

Spectrosc. 238 (2006) 23-40. 
[15] L. Haloncn, T. Carrington Jr, J. Chcm. Phys. 88 (1988) 4171- 

4185. 

[16] F. Iachello, S. Oss, J. Mol. Spcctrosc. 142 (1990) 85-107. 
[17] Y. Zheng, S. Ding, J. Mol. Spectrosc. 201 (2000) 109-115. 
[18] O. Naumenko, A. Campargue, J. Mol. Spectrosc. 210 (2001) 
224-252. 

[19] Vl.G. Tyutercv, Chcm. Phys. Lett. 348 (2001) 223-234. 

[20] Vl.G. Tyuterev, L. Rcgalia-Jarlot, D.W. Schwenkc, S.A. 

Tashkun, Y.G. Borkov, C. R. Phys. 5 (2004)189-199. 
[21] F. Michelot, J. Moret-Bailly, J. Physique 48 (1987) 51-72. 
[22] C. Leroy, F. Michelot, J. Mol. Spcctrosc. 151 (1992) 71-96. 
[23] C. Leroy, F. Michelot, Can. J. Phys. 72 (1994) 274-289. 
[24] F. Michelot, M. Rcy, Eur. Phys. J. D, 30 (2004) 181-189. 
[25] F. Michelot, M. Rey, Eur. Phys. J. D, 33 (2005) 357-386. 
[26] F. Michelot, M. Rcy, Eur. Phys. J. D, 44 (2007) 467-495. 
[27] B.C. Wybournc, Classical Lie Groups for Physicists, Wiley- 

Intcrscicnce, New York, 1974. 
[28] I.M. Gcl'fand, M.L. Zetlin, Dokl. Akad. Nauk. (in Russian) 71 

(1950) 825-828. 
[29] J. Louck, Amer. J. Phys., 38 (1970) 3-42. 

[30] A.R. Edmonds, Angular Momentum in Quantum Mechanics, 

Princeton University Press, 1974. 
[31] W. Holmann III, L.C. Bicdenharn, in Group Theory and Its 

Applications, edited by E.M. Loebl, Academic Press, New York, 

1971 Vol. II, p. 1-73. 
[32] C. Leroy, PhD thesis, University of Dijon, France 1991. 
[33] V. Boujut, PhD thesis, University of Dijon, France 1996. 



[34] L. Pluchart, C. Leroy, N. Sanzharov, F. Michelot, E.S. 
Bekhtereva, O.N. Ulenikov, J. Mol. Spectrosc. 232 (2005) 119- 
136. 

[35] N.A. Sanzharov, C. Leroy, O.N. Ulenikov, E.S. Bekhtereva, J. 

Mol. Spcctrosc, 247 (2008) 1-24. 
[36] L.R. Peebles, P. Marshall, J. Chcm. Phys. 117 (2002) 3132-3138. 



Tabic 3: Observed and calculated energies of D2S with 22 experimental data. (1 < P < 12). 
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Tabic 3: (cont.) 
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Tabic 3: (cont.) 
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Tabic 3: (cont.) 
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Tabic 3: (cont.) 
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